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. Abstract 

O 

I We present analytic expressions to evaluate at 0{a) the effects of soft-photon emission, 



IQ ■ and the related virtual corrections, in non-leptonic decays of the type B,D ^ P1P2, 
where Pi^2 are scalar or pseudoscalar particles. The phenomenological implications of 
these results are briefly discussed. For B decays into charged pions the effects of soft- 
I , photon emission are quite large: the corrections to the rates can easily exceed the 5% 
level if tight cuts on the photon energy are applied. 



O 



X : 1 Introduction 



The large amount of data collected at B factories has allowed to reach statistical accuracies 
of the order of a few percent on the measurements of several B meson branching fractions. 
At this level of accuracy electromagnetic effects cannot be neglected. On the one hand, in 
order to ensure a good control of the experimental efficiencies, it is necessary to include 
in Montecarlo simulations the unavoidable emission of soft photons that accompanies all 
processes with charged particles. On the other hand, the effective cuts applied on the 
(soft) photon spectra are a key information for a meaningful comparison between theory 
and experiments. 

The theoretical treatment of the infrared singularities generated within QED is a 
well known subject and one of the pillars of quantum field theory. A clear and very 
general discussion can be found, for instance, in the classical papers |2] . The general 
properties of QED have been exploited in great detail in the case of genuine electroweak 
processes, or processes which can be fully described within perturbation theory within the 
Standard Model (SM). In these cases there exist both precise theoretical calculations of 
the electromagnetic (e.m.) corrections and accurate Montecarlo programs which include 
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the effects of soft photon emission, such as PHOTOS [B^. Similar tools have not been 
developed for most exclusive hadronic processes and, in particular, for B and D decays. 

Recently, the issue of electromagnetic corrections have received considerable attention 
in the context of kaon physics [H El El EI- As discussed in these works, and as confirmed 
by recent experimental analyses |B] , a correct simulation of electromagnetic corrections is 
a key ingredient for a precise determination of Vus and other effective couplings of weak 
interactions. 

The purpose of the present work is to present simple analytic formulae for he theo- 
retical evaluation, and the numerical simulation, of the leading radiative corrections in B 
or D meson decays into two scalar or pseudoscalar particles. Given the universal char- 
acter of the infrared singularities, we evaluate the effects of soft photon emission (and 
the corresponding virtual corrections) within scalar QED and in the approximation of a 
point-like effective weak vertex. The results thus obtained are valid up to constant 0{a) 
terms (not enhanced by large logs) related to the matching between this effective theory 
and the "true" theory where the dynamical aspects of weak interactions are taken into 
account. Moreover, our calculation does not take into account the possible 0{E^) terms 
associated to non-bremsstrahlung amplitudes (hard photon emission). 

2 Photon-inclusive widths 

The most convenient infrared-safe observable related to the process H P1P2 is the 
photon-inclusive width 

r-i(E— ) = T{H^ P,P, + n^) I J. , (1) 

namely the width for the process H P1P2 accompanied by any number of (undetected) 
photons, with total missing energy less or equal to E'^^^ in the H meson rest frame. At 
any order in perturbation theory we can decompose T^{^^ in terms of two theoretical quan- 
tities: the so-called non-radiative width, r^2^ the corresponding energy-dependent 
e.m. correction factor Gi2{E^^^), 

pincl^^max^ = T% G^iE'^'"") ■ (2) 

The energy dependence of GuiE) is unambiguous and universal (i.e. independent from the 
short-distance dynamics which originate the decay) up to terms which vanish in the limit 
E ^ 0. On the contrary, the normalization of G 12(E) is arbitrary: we can always move 
part of the finite (energy- independent) electromagnetic corrections from to Gi2{E). 
Only the product in corresponds to an observable quantity. 

In the following we report the explicit expressions of G 12(E) as obtained by means of a 
scalar-QED calculation at 0{a). In order to treat separately infrared (IR) and ultraviolet 
(UV) divergences, we regulate the former by means of a photon mass and the latter by 
means of dimensional regularization. We then renormalize the point-like weak coupling 
in the MS scheme. 
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By construction, we define the non-radiative amplitude as follows 
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namely the tree-level rate expressed in terms of the renormalized weak coupling. With 
this convention, the function Gi2{E) can be written as 
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where, following the notation of Ref. j^, we have denoted by H12 the finite term arising 
from virtual corrections, and by F12 the energy-independent term generated by the real 
emission: 
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As expected, after summing real and virtual corrections the infrared logarithmic diver- 
gence cancel out in Gi2{E), giving rise to the universal hi{MH/E) terms proportional 
to 
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Note that Gi2{E) does depend explicitly on the ultraviolet renormalization scale /x. 
The scale dependence contained in Nui^fi) cancels out only in the product r'j'2 Gui^E) due 
to the corresponding scale dependence of the weak coupling. As we shall discuss in more 
detail later on, in practice we are not able to exploit the numerical consequences of this 
cancellation due to the absence of a first-principle calculation of AH^PiP2{f^)- 

In the — > P1P2 case we find the following explicit expressions for the coefficients 
in eq. (0): 
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where 
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In the if^ — > P^^ P2 case we find: 
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In both cases {H^ Pi P2 ^^d if ^ Pj^Pj^) fully recover the results of Ref. |1] in 
the limiting case ri = r2 and setting /i = Mi. 



3 Numerical results and differential distributions 

In figure Q we report a few examples of the e.m. correction factors relevant to P — >■ 
TTTT, TiK, and KK modes. As can be noted, the effects are quite large in case of tight 
cuts on the photon energy or, equivalently, tight cuts on the invariant mass of the two 
pseudoscalar mesons. The e.m. effects grow logarithmically with respect to the velocities 
of the charged particles in the final state and with respect to the cut on the maximal 
photon energy (relative to the mass of the decaying particle). As a result, the largest 
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Figure 1: Examples of the e.m. correction factors for the photon- inclusive widths as a 
function of the cut on the maximal missing energy. Left: Gtt+tt- (full lines) and Gk+k- 
(dashed lines) for the corresponding B decays; in each set the two curves are obtained for 
11 = Mp (upper curve) and /i = M^(Mi^) (lower curve). Right: scale-independent ratios 
Gko-k+ /Gk+ttO (full line) and Gt^+t,- /Gj^+t^o (dashed line) for the corresponding B decays. 



effects are found for B decays into charged pions. For instance, the observable ratio 
r[i? — > 7r^7r^(7)]/r[i? — > tt^tt^], determined with a cut of about 250 MeV on the missing 
mass of the tt+tt" pair, receives a 5% negative isospin-breaking correction with respect to 
the corresponding theoretical ratio evaluated in absence of soft-photon contributions. 

As anticipated, the overall normalization of the GulE) factors is convention depen- 
dent. By construction, they include the whole effect of the bremsstrahlung but only part 
of the virtual e.m. corrections, namely the leading soft contribution which can be reliably 
estimated within scalar QED. The separation between hard and soft virtual contributions 
is arbitrary and determined by the choice of the UV renormalization scale /x. The most 
natural choice for the latter is between the light pseudoscalar masses (M^r or M^) and 
Mp ^ 770 MeV. As shown in the left plot of figure [H the uncertainty associated with this 
range is at most of (9(1%) and definitely subleading with respect to the leading effects 
induced by soft contributions. 

In principle, this scale uncertainty could be removed by means of an appropriate 
matching between our effective theory -where photons and mesons are the only dynamical 
degrees of freedom and each weak amplitude is determined by a momentum-independent 
coupling- and more predictive effective theories where the weak decay amplitudes are 
computed in terms of the fundamental SM couplings. In the last few years there has been 
a substantial progress toward a precise estimate of non-leptonic weak decay amplitudes 
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within the SM (see e.g. Ref. jHj and references therein); however, we are still far from the 
(9(1%) level, especially for the decays into two light pseudoscalar mesons. Moreover, all 
present approaches do not include the structure-dependent electromagnetic corrections 
which would allow to match the scale dependence of the G 12(E) factors. 

Despite the absence of a precise UV matching, in a few cases the e.m. correction factors 
we have computed allows to evaluate in a precise way the amount of isospin breaking 
induced by soft photons. As illustrated in the right plot of figure ^ the scale uncertainty 
drops out completely in all the ratios of two-body photon-inclusive widths with at least 
one charged particle in the final state. 

It is worth to stress that the e.m. corrections discussed in this work can be one of the 
sources of the "anomalous" isospin-breaking effects identified in recent phenomenological 
analyses of B ^ ttK and B ^ tttt decays PU] • A consistent phenomenological analysis of 
these channels with the inclusion of radiative corrections requires the information on the 
experimental cuts applied on the soft-photon radiation and is beyond the purpose of the 
present work. However, it is clear that radiative corrections need to be included in these 
channels, especially if one is interested in isospin-breaking effects as a tool to identify 
possible deviations from the SM. 

From the purely experimental side, the only relevant aspect of the Gi2{E) factors is 
their energy dependence, which is unambiguously determined up to 0{E) terms. This 
allows us to evaluate the observable missing-energy distribution, or the soft-photon spec- 
trum, dr^^2^{E)/dE. As is well known, the E —>■ singularity of this distribution is not 
integrable if evaluated at any fixed order in perturbation theory; however, the all-order 
resummation of the leading infrared singularities leads to an integrable distribution 
In our case, this can be written as 
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with the hi2 coefficients given in eq. ((Tj). The exponentiation of the leading- log corrections 
does not lead to appreciable numerical difference with respect to the pure 0{a) result 
for missing energies above few MeV. However, as discussed in [7j, having an integrable 
distribution can be very useful in preparing a dedicated Montecarlo program. 

Concerning the angular distribution of the bremsstrahlung photons, the 0{a) result 
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where E^ and 9^ denote, respectively, photon energy and angle between photon and Pi 
momenta in the if meson rest frame, and 
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The Ri2 coefficients, defined by 
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where Qi^2 and ~Qh denote, respectively, the electric charges of Pi^2 and H in units of 
the electron charge, assume the following explicit form 
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4 Conclusions 

In the last few years there has been a substantial progress in the experimental determi- 
nation of two-body non-leptonic B decays into light pseudoscalar mesons. Several new 
theoretical tools have also been developed for the description of these interesting pro- 
cesses. However, in most cases electromagnetic effects of long-distance origin have been 
ignored, both in the theoretical predictions and in the experimental analyses. Aiming at 
accuracies of the order of a few percent, this approximation is no longer valid. 

In this work we have presented a detailed discussion of the electromagnetic effects of 
long-distance origin in two-body non-leptonic B decays. In particular, we have computed 
the leading 0{a) effects induced by both real and virtual photons in a generic process of 
the type H — > PiP2{'y), where both H and Pi^2 are scalar or pseudoscalar particles (with 
arbitrary masses). The results thus obtained can be applied to both B and D decays. 

We have performed the calculation within scalar QED and assuming an effective con- 
stant coupling for the weak vertex. Given the universal character of the infrared singu- 
larities, this simplified framework allows to identify the leading effects induced by soft 
photons. The results obtained within this framework are valid up to energy-independent 
0{a) terms related to the matching between this effective theory and the SM, and possi- 
ble 0{E^) terms associated to the emission of hard photons (which should be identified 
experimentally). In the case of the photon-inclusive rates, the analytic formulae we have 
obtained represent a generalization of the results of Ref . Q , which are recovered in the lim- 
iting case Ml = M2. In order to facilitate the experimental simulation of the soft-photon 
radiation, we have also discussed the corresponding angular and energy distribution. 

From the phenomenological point of view, the largest impact of the soft-photon radi- 
ation is found in the B tt+tt" channel, where the corrections to the partially-inclusive 
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rate can easily exceed the 5% level. We stress that detailed estimate of these effects can 
only be made once the experimental information on the photon cuts (or the missing en- 
ergy) is available. Without this information, an unambiguous comparison between theory 
and experiments, and also the combination of different experimental results, cannot be 
performed. 
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